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On Tightly Bounding the Dubins Traveling Salesman’s Optimum 
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Abstract 

The Dubins Traveling Salesman Problem (DTSP) has generated significant interest over the last decade 
due to its occurrence in several civil and military surveillance applications. Currently, there is no algorithm 
that can find an optimal solution to the problem. In addition, relaxing the motion constraints and solving 
the resulting Euclidean TSP (ETSP) provides the only lower bound available for the problem. Plow ever, in 
many problem instances, the lower bound computed by solving the ETSP is far below the cost of the feasible 
solutions obtained by some well-known algorithms for the DTSP. This article addresses this fundamental 
issue and presents the first systematic procedure for developing tight lower bounds for the DTSP. 


1. Introduction 

Given a set of targets on a plane and a constant p > 0, the Dubins Traveling Salesman Problem (DTSP) 
aims to find a path such that each target is visited at least once, the radius of curvature of any point in the path 
is at least equal to p, and the length of the path is minimal. This problem is a generalization of the Euclidean 
TSP (ETSP) and is NP-hard |[^[^. The DTSP belongs to a class of task allocation and path planning 
problems envisioned for a team of unmanned aerial vehicles in Q. The DTSP has received significant 
attention in the literature ||T|[2]|4]-[T^, mainly due to its importance in unmanned vehicle applications, the 
simplicity of the problem statement, and its status as a hard problem to solve because it inherits features 
from both optimal control and combinatorial optimization. 

Currently, there is no procedure for finding an optimal solution for the DTSP. Therefore, heuristics and 
approximation algorithms have been developed over the last decade to find feasible solutions. Tang and 
Ozguner | [TT[ present gradient-based heuristics for both single and multiple vehicle variants of the DTSP. 
Savla et al. use an optimal solution to the ETSP to find a feasible solution for the DTSP, and they bound 
the cost of the feasible solution with respect to the optimal cost of the ETSP. Rathinam et al. Q develop 
an approximation algorithm for the DTSP in cases where the distance between any two targets is at least 
equal to 2p. Ny et al. Q develop an approximation algorithm for the DTSP in which the approximation 
guarantee is inversely proportional to the minimum distance between any two targets. The weakness of the 
approximation guarantees of these algorithms for the DTSP is due to the lack of a good lower bound, as all 
these algorithms essentially use the Euclidean distances between the targets to bound the cost of a feasible 
solution. 

Other heuristics have been used for solving the DTSP. A receding horizon approach that involves finding 
an optimal Dubins path through three consecutive targets is used to generate feasible solutions in p4| . The 
heuristic in Q finds a feasible solution by minimizing the sum of the distances travelled by the vehicle and 
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the sum of the changes in the heading angles at each of the targets. Macharet et al. l|^|^ first obtain a tour by 
solving the ETSP and then select the heading angle at each target using an orientation-assignment heuristic. 
A multiple lookahead approach is used to find feasible solufions in |[^[T^. Mefa-heurisfics have also been 
developed fo find feasible solufions for fhe DTSP in l|9 101. 



Figure 1. There are four possible headings at each target. A feasible solution for the DTSP can be 
obtained by choosing a heading at each target and finding a corresponding optimal TSP path. 


Another common approach pA][T^ involves discretizing the heading angle at each target and posing the 
resulting problem as a one-in-a-set TSP (Fig. [^. The greater the number of discretizations, the closer an 
optimal one-in-a-set TSP solution gets to the optimal DTSP solution. This approach provides a natural way 
to find a good, feasible solution to the problem p6| . However, this also requires us to solve a large one-in- 


a-set TSP, which is combinatorially hard. Nevertheless, this approach provides an upper bound 1131 for the 
optimal cost of the DTSP, and simulation results indicate that the cost of the solutions start to converge with 
more than 15 discretizations at each target. 

The fundamental question with regard to all the above heuristics and approximation algorithms is how 
close a feasible solution actually is to the optimum. For example. Fig. shows the cost of the feasible 
solutions obtained by solving the one-in-a-set TSP and the ETSP for 25 instances with 20 targets in each 
instance. Even with 32 discretizations of the possible angles at each target, the cost of the feasible solution is 
at least 30% greater than the corresponding optimal ETSP cost for several of these instances. As the optimal 
cost is not known for the DTSP, identifying a tight lower bound is crucial for determining the quality of the 
solutions that have been provided as well as for developing constant factor approximation algorithms. 

In these algo- 


This fundamental question was the motivation for the bounding algorithms in 117 


rithms, the requirement that the arrival and departure angles must be equal at each target is removed, and 
instead there is a penalty in the objective function whenever the requirement is violated. This results in a 
max-min problem where the minimization problem is an asymmetric TSP (ATSP) and the cost of travel¬ 
ing between any two targets requires solving a new optimal control problem. In terms of lower bounding, 
the difficulty with this approach is that we are not currently aware of any algorithm that will guarantee a 
lower bound for the optimal control problem. Nonetheless, this is a useful approach, and advances in lower 
bounding optimal control problems will lead to finding lower bounds for the DTSP. 

In this article, we propose a new approach to finding tight lower bounds for the DTSP. This is the first 
systematic procedure available for the DTSP and is a natural counterpart to the one-in-a-set TSP approach 
we discussed above. In this approach, we remove the requirement that the arrival angle and the departure 
angle at each target must be the same, but we restrain these angles so that they belong to one sector or 
interval (refer to Fig. |^. The lower bounding problem aims to choose an interval at each target such that 
the arrival angle and the departure angle at the target belong to the same interval, each target is visited at 









Instance Number 


Figure 2. A comparison between the cost of the feasibie soiution (upper bound) obtained by soiving 
the one-in-a-set TSP with 32 discretizations and the optimai cost of the corresponding Euciidean 
TSP (iower bound) for 25 instances. There are 20 targets in each instance, and the iocation of each 
target is sampied from a 1000x1000 square. Aiso, the minimum turning radius of the vehicie is set 
to 100. 


least once, and the sum of the costs of travelling between the targets is minimized. The cost of traveling 
between two intervals corresponding to two distinct targets now reduces to a new optimal control problem, 
which we refer to as the Dubins interval problem. Given two targets and an interval at each target, the 
problem is to find a Dubins path such that the departure angle at the initial target and the arrival angle at 
the final largel belong fo fhe given infervals and fhe lengfh of fhe pafh is minimal. The lower bounding 
problem is a one-in-a-sef TSP and can be solved jusf like fhe upper bounding problem. If fhe size of each 
of fhe infervals af each largel reduces lo zero, fhe lower bounding problem reduces fo fhe DTSP. If Ihere is 
only one interval of size In al each largel, fhe resull is a Euclidean TSP. As fhe size of fhe infervals al fhe 
largels becomes smaller, fhe one-in-a-sef TSP becomes combinaforially hard, similar lo fhe upper bounding 
problem. Neverlheless, Ihis provides a sysfemafic approach fo finding lower bounds for fhe DTSP, provided 
fhe Dubins inlerval problem can be solved. 

The Dubins interval problem is a new generalizafion of fhe slandard Dubins problem |25| which has 
nof been formulaled or solvecQin fhe liferalure. The difficully wilh solving Ihis problem lies in fhe facf fhaf 
the length of the shortest Dubins paths between any two targets is a non-linear, discontinues function of 
the heading angles of the targets. Therefore, finding fhe optimal heading angles from fhe given infervals al 
fhe largels lhal minimizes fhe lengfh of fhe Dubins pafh is non-lrivial. In Ibis article, we solve fhe Dubins 
inlerval problem using fhe monolonicily properlies and fhe exlremal values of fhe lengfh of fhe Dubins palhs. 


The following are fhe confribulions of Ihis article: 

*We point out that we are also working on a completely different approach which is based on optimal control in j^. However, 
the approach in (D still relies on the results of this paper. Second, unlike the analysis and results in this paper, the work in (ID does 
not provide information about the rate of change of lengths of the Dubins paths as a function of the heading angles at the targets. 
This information is critical and very useful in the development of bounds for feasible solutions and approximation algorithms for 
the DTSP. 






Figure 3. There are four intervals at each target. A lower bound for the DTSP can be obtained by 
choosing an interval and restricting both the arrival and the departure angles to be in the chosen 
interval at each target, and then finding a corresponding optimal TSP path. The shaded interval at 
each target shows the chosen interval with the arrival and departure angles in blue. 


1. The formulation of the lower bounding problem for the DTSP as a novel one-in-a-set TSP where the 
cost of traveling between any two targets requires solving a Dubins Interval Problem. This is the first 
formulation that aims to provide a tight lower bound for the DTSP 

2. The first algorithm to solve the Dubins Interval Problem by exploiting its structure and monotonicity 
properties. 

3. Numerical results to corroborate the performance of the proposed lower bounding approach for the 
25 instances shown in figure]^ 

2. Lower Bounding Problem Formulation 


The set of targets is denoted by T = {1,2, • • • ,«}, where n is the number of targets. The set of avail¬ 
able angles [0,27r] at any target i is partitioned into a collection of closed intervals denoted by := 
{[ 0 ,(Pii],[(Pii,(Pi 2 ],-■ ■ , )<Pmi, = 271]}, whcrc mi{> 1) denotes the number of intervals at target i and 

the (pij are constants such that 0 < (pn < (pa < • < (Piim = 271. Let {xi,yi) denote the location of target 

i G T, and let the arrival angle and the departure angle of the vehicle at target i be denoted by and 
respectively. The configuration of the vehicle leaving target i at is then denoted by (x,,y/, 6,y), and 
(x/,y;, dia) similarly denotes the vehicle’s arrival configuration. The length of the shortest Dubins path from 
(x/,y/, did) to {xj,yj, Oja) is denoted by dij{6id, Oja)- Given an interval /,■ at target i and an interval Ij at target 
j, define d*j{Ii,Ij) := dij (^Oid^ Oja)- The objective of the Bounding Problem (BP) is to find a 

sequence of targets (^i ,^ 2 , • • • ,i'n), i'/ G T, to visit and an interval Is- G for each target si G T such that 


• each target is visited at least once, and 

• the cost I”-/ d*s.s.^^ {Is,, ) + dt„si (2i„, 4,) is minimized. 

Addressing this BP first requires solving £.^^i.dij{0id,0ja)- Once this problem is solved, the 

BP is essentially a one-in-a-set TSP. In this article, we transform the one-in-a-set TSP into an ATSP using 


the Noon-Bean transformation |22| and then convert the resulting ATSP into a symmetric TSP using the 


transformation in |23|. The symmetric TSP is solved using the Concorde solver |24| to find an optimal 
solution. 

Prior to addressing the Dubins Interval Problem in the next section, we first formally state the lower 
bounding result in the following proposition. 






Proposition 2.1, The optimal cost to the BP is a lower bound to the DTSP 

Proof. Any optimal solution to the DTSP is a feasible solution to the BP for any positive number of intervals 
at each target. Therefore, the optimal cost of the BP must be a lower bound to the optimal cost of the 
DTSP. □ 


3. Dubins Interval Problem 


Without loss of generality, let the Dubins interval problem be denoted as minej^/j where 

d\2{Q\,Q2) indicates the shortest path (also referred to as the Dubins path) for traveling from (xi,yi,6i) 
to (x2,y2,S2) subject to the minimum turning radius constraint (Fig. |^. Here the interval 4 is defined 
as C [ 0 , 27 r] for k = 1 , 2 . Given an initial configuration (xi,yi,6i) and a final configuration 

(x2,y2, 02)> L.E. Dubins |251 showed that the shortest path for a vehicle to travel between the two configu¬ 
rations subject to the minimum turning radius (p) constraint must consist of at most three segments, where 
each segment is a circle of radius p or a straight line. Specifically, if a curved segment of radius p along 
which the vehicle travels in a counterclockwise (clockwise) rotational motion is denoted by L(R), and the 
segment along which the vehicle travels straight is denoted by S, then the shortest path is one of RSP, RSL, 
LSR, LSI, RLR, and LRL. 



Figure 4. A feasible solution to the Dubins interval problem. 

Let RSL{d \, 62) denote the length of the RSL path from {x\,y \, 0 i) to (x2,y2, ^2)- RSL{di , 62) is set to 00 
if the RSL path does not exist. Let RSR{6i , 62), LSR{6i , 62), LSL{d \, 82), RLR{d \, 82), and LRL{8\ ,82) be 
defined in a similar way. Using these definitions, the Dubins interval problem can be written as follows: 


min di2{8i,82)= min {RSR(8u82),RSL(8i,82),LSR(8u82),LSL(8i,82),RLR(8i,82),LRL(8i,82)}. 

016/1,02 eh 016/1,02 6/2 

( 1 ) 

Remark 3.1. <ii2(0ij02) A a lower semicontinuous function and is minimized over closed and bounded 
intervals 1 \ and 4. Therefore, the Dubins interval problem is well defined, i.e., there exist 8^ G 4 and 
82 G 4 such that di2{8*, 82) = ‘^12(81,82). 

To solve the Dubins interval problem, we also consider shortest paths that contain at most two segments 
between {x\,y\) and (x2,y2)- For any path E {RS,LS,SR,SL,RL,LR} and 81 E 4 , let ^^{81) denote 
the distance of the shortest path of type ^ that starts at (xi,yi) with a departure angle of 61 and arrives at 



{x2,y2) with an arrival angle in I2. In this case, the arrival angle at (^2,^2) will be a function of 61 and .!7 and 
is denoted as 62 (=^ 7 ,61)■ =^^(61) is set to 00 if a path of type ^ does not exist or if 62 (=^ 7 , ^i) ^h- Similarly, 
let ^^{62) denote the distance of the shortest path of type ^ that starts at (xi,yi) with a departure angle 
in I\ and arrives at {x2,y2) with an arrival angle of 62- In this case, the departure angle at (xi,yi) will be a 
function of 62 and ^ and is denoted as 61 (^, 02)- ^^{^2) is set to 00 if the path of type .!7 does not exist 
or if 01 , 62) ^ /i- From the definitions, note that mingi^/j ^^(di) = mine, ^^^(62). 

The following theorem provides a way to further simplify equation ([T]) and solve the Dubins interval 
problem: 

Theorem 3.1. 

min min{(fi 2 ( 0 i, 02 )} = min{(jf*,min{/?5'^(0i),5'/?V0i),L5'^(0i),5'L'(0i),L/?V0i),/?LV0i)}} (2) 

016/1 e/2 016/1 


where 


d* :=mm{dn{ef", 


02"™),/fi2(0r“,02™”),//i2(er 


02™“),//i2(0r“, 


02"'“)}. 


In English, this theorem states that an optimal path to the Dubins interval problem must be one of the 
following: 

1 . An optimal Dubins path consisting of at most three segments such that both the arrival and departure 
angles at each target belong to one of the boundary values of the respective intervals, or 

2 . An optimal Dubins path consisting of at most two segments such that the angle constraints are satis¬ 
fied. 


After proving this theorem in the next subsection, we will provide algorithms to solve for the optimal 
Dubins paths with at most two segments (i.e., to solve min^ig/j ^( 0 i) for any path ^ G {RS^ ,SR^ ,LS^ ,SL^ ,LR^ ,RL^}). 
As d* in the above theorem can already be computed directly using Dubins’s result p 5 | , one can then com¬ 
pute the optimal cost for the Dubins interval problem and the corresponding departure and arrival angles. 

3.1. Proof of Theorem |3.1 1 

This theorem will be proved in two parts. First, we will first show how to simplify min^, 6/1,026/2 >^(^11^2) 
for any path G {RSR,RSL,LSR,LSL}. Then we will address the LRL and the RLR paths. 


3 . 1 . 1 . Optimizing RSR, RSL, LSR, and LSL paths. The following result is known [261 for each of the paths 
G {RSR,RSL,LSR,LSL} from (xi,yi, 0 i) to (x2,y2,02): 


Lemma 3.1. For any ^ G {RSR,RSL,LSR,LSL} and i= 1,2, either > 0 V 0; or < 0 

V 0 ( when exists and none of its curved segments vanish. 


Now let us apply the above lemma to the RSL path. The RSL path ceases to exist when the segment S 
vanishes, i.e., the RSL path reduces to an RL path. In addition, when one of the curved segments vanishes, 
the RSL path reduces to either the RS or the SL path (refer to Fig. |^. Therefore, given d\, the optimum for 
ming^ 6 [ 0 ”"' 01 '“] RSL{ 6 \, 62 ) must be attained when 62 = 6 f“" or 02 = 0 ™“ or when the RSL path reduces to 
an RL, RS, or SL path. This can be stated as follows: 
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Figure 5. Given Ou the iength of the RSL path varies monotonicaiiy with respect to O 2 wherever the 
path exists and none of its curved segments vanish. 


min{/?5L(ei,02)} := mm{RSL{euef'^),RSL{eue^"^),RS\ei),SL\ei),RL\Gi)}. (3) 

02 e/2 


Therefore, 


min min{/?5'L(0i,02)} 

0i6/i02e/2 

= min mm{RSL{di, df"),RSL{di, 6^“^),RS\ei),SL\ei),RL^ {di)} 

0ie/i 

= mini min RSUdi, dT "), min RSUdi , 02™“), min {RS^ (01), (0i), (0i)}}. (4) 

016/i 016/i 01 e/i 

Similarly, using Lemma [TT] again, we get the following: 

mmRSL{euer)=mm{RSL{d'r,ef'^),RSL{er\eT‘"),RS^(.Q2n,SL\dr),RL^{ef")}, (5) 

0ie/i 

min/?SL(0i, er"") = min{7?SL(er”, 02™“),7?SL(e;”“, e^^^),RS^{e^‘^),SL^{eT‘^’‘),RL^{e^‘^)}. (6) 

0l6/i 


Now, one can easily verify the following: 


For any £r e {RS,SL,RL}, min and min (7) 

0l6/i 0i6/i 


Substituting for min^j g/j RSL{di , 0™'”) and min^jg/j RSL{di , 0™^) in Q using equations Q and Q and 
simplifying further using Q, we get 


min min7?SL(0i, 02) = mm{RSL*, min{/?S^(0i),SL^(0i),/?L^(0i)}}, 

Oi^Ii 92^12 9i^Ii 


( 8 ) 



where 


RSL* := mm{RSL{df’^, ef”),RSL{e'^‘^\ ef'^),RSL[e'r, 6!^“’^),RSL[e'r\ 02"'“)}- 

As Lemma[37T]is also applicable to RSR, LSL, and LSR paths, one can use the above procedure and sim¬ 
plify min0;g/j_e2g/2 7?S/?(ei,02), min0;e/j_e2g/2LSL(ei,62), and min0je/j_02e/2L5'7?(ei,62) in a similar way. 
Combining all these results, we obtain the following: 


min min min ^(61,62) = min{,^*, min{/?5'^(0i),5'/?'(0i),L5'^(0i),5'L'(0i),L/?'(0i),/?L^(0i)}} 

di&li ^e{RSR,RSL,LSR,LSL} 0l6/i 

(9) 

where 


; = 


min mm{^(er'', 

3^€{RSR,RSL,LSR,LSL} 


02™’”), ^(0r“, 02™'”), ^(0;”'”, 02"““), ^(0r'“, 


02"’“) }. 


3.1.2. Optimizing RLR and LRL paths. Xavier et al. 1261 have shown that the RLR and LRL paths can¬ 
not lead to an optimal Dubins path if the distance between the two targets is greater than 4p. Therefore, 
in this section, we assume that the distance between the two targets is at most 4p. We will focus on 
min0jg/j 02g/2L/?L(0i,02); min0,e/j 02e/2/?L/?(0i,02) can be solved in a similar way. Given 0i, unlike the 
length of the RSL path, LRL{6i ,62) is not monotonous with respect to 62 when LRL exists. Without loss of 
generality, we assume that 0i = 0 and first aim to understand LRL{0,62) as a function of 62 (refer to Fig. [^. 
Target 1 is located at the origin and target 2 is located at (x,y). The angles a and j3 in Fig. |^are functions 
of 02. For brevity, we use a and j8 in place of a(02) and j3(02), respectively. Let LRL(0,62) be denoted as 
®(02) := (27r-|-2a-|-2j8 -|- 02)p. In the ensuing discussion, we use the fact that the length of the R segment 
in an LRL path must be greater than Tip (i.e., 0 < a -|- j8 < tt) for the LRL path to be an optimal Dubins path 
between any two targets |25 27|. 



Figure 6. LRL path for 9 i = 0. 

Lemma 3.2. If the LRL path exists and none of its curved segments vanish, then for any 02 such that 
0 < a{d2) + ji (62) < 71 , 1^/0 except when T)(02) reaches a maximum, when 62 satisfies a + 7 t /2 = 62- 









Proof. Using Fig. a and j 3 can be obtained in terms of 62 as follows: 


2 psina + p = 2 psinj 3 +pcos02+Jj (10) 

2p cos a + 2p cos j8 + p sin 62 =( 11 ) 

Differentiating and simplifying the above equations, we get 


da „ dB sm62 

dd2 ^ dd2 2 

da . n cos 62 

sma—+sm/i—= —, 


Further solving for the derivatives, we get 


d [5 cos(62 — a) 

dd2 2 sin(a + j 8 )’ 

da cos(02 + i 3 ) 

dd2 2 sin(a + j 8 )’ 


Therefore, 


dD 

dd2 


,^dB ^da 


p( 


COS(02 — of) 


I COs{d2+P) ^ 


sin(a + j 8 ) sin(a + j 3 ) 


( 12 ) 

( 13 ) 

( 14 ) 

( 15 ) 

( 16 ) 
( 17 ) 


Equation = 0 yields the following possibilities: 62 = f + o;or62 + i8 = —f. 62 + P = —j corresponds 
to the case where the second left turn disappears; there is a jump in the length of the LRL path at this 62, 
and therefore does not exist. 62 = j + a corresponds to the case where the turn angle in the right turn 
is equal to the turn angle in the second left turn; one can verify that ^(62) reaches a maximum at this point 
because ^ < 0 (refer to Fig. B. 

□ 


The derivatives of L 7 ?L( 6 i,62) do not exist when any turn in the path disappears or when the angle in 
the right turn becomes equal to tt, as shown in Fig. The length of the two paths (Fig. when the LRL 
path just ceases to exist are denoted by L/?L^( 0 i) and L/?L^( 0 i). Therefore, applying the above lemma to 
the LRL path and following similar steps to those in subsection | 3 . 1 . 1 [ we get the following result: 


min{L/?L( 0 i,e 2 )} := min{LRL(di,df'^),LRL(ei,dr^),LR\di),RL\di),LRL^(di),LRLl(di)}. ( 18 ) 

02 6/2 

Again, as in subsection | 3 . 1.1 [ one can further simplify the above optimization problem: 

min min {LRLfO,,62)} = min{LRL*, mm{LR\di),RL\di),LRL\(df),LRLI ()}} ( 19 ) 

01 6/1 02 6/2 01 6/1 


LRL* := mm{LRL{eff df''),LRL{e’f^, "),L/?L(ej” 


, ef^),LRL{e’f 


Qmax^Y 


where 














82 


Figure 7. Given 0i, the iength of the LRL path reaches a maximum when 02 = f + «, as shown. This 
figure aiso shows the vaiues of 02 where the LRL path just ceases to exist. 


Note that LRL\{Q\) and LRL\{d\) can never result in an optimal Dubins path because the angle in the 
right turn is equal to n p7| . Therefore, once equation ( [T^ is substituted in equation ([T]l, the functions 
LRL\{di) and LRL\{d\) will drop out. 

mine|g/j mine2g/2{/?L/?(0i, 02)} can be simplified in a similar way. Hence, combining the above results 
with equation Q, we obtain the result stated in Theorem 

4. Algorithms for Optimizing Dubins Paths with At Most Two Segments 

The only remaining step needed to solve the Dubins interval problem is to show how to optimize 
min0,e/j{7?S^(0i),57?'(0i),L5'^(0i),SL'(0i), LR^{Q\),RL}{Q\)}. In this section, we will solve two prob¬ 
lems: min0jg/|{7?S^(0i)} and mine, g7j{/?L^(0i)}. The remaining paths can be solved using simple transfor¬ 
mations (reflections about the a: or the y axis). 

4.1. Optimizing the RS path 

Without loss of generality, a reference frame can be chosen such that target 1 is at the origin and target 2 
lies on the x-axis as shown in Fig. Here x represents the Euclidean distance between the targets. Given 0i, 
the existence of the RS path as well as its length can be determined using geometry. The length of the S path, 
the angle between the x-axis and the S path, and the final arrival angle at target 2 are also functions of 0i 
and can be expressed as L{d\), 0(0i), and 62 {RS, 0i), respectively. Let the length of the RS path be denoted 
as T)(0i). For brevity, in some places we will use L, 0,02, and D instead of L(0i),0(0i),02(/?5', 0i), and 
T)(0i), respectively. Let ds : = x if the angle of the straight line joining the two targets lies in the intervals 
I\ and h- If the angle constraints are not satisfied, ds is set to oo. Similarly, let d^ denote the length of the 
shortest circular arc of type R that joins the two targets such that the boundary angles of the arc belong to 
the respective intervals at the targets. If such an arc does not exist, dR is set to oo. In the following lemma, 
we assume that [0™",0f“] C [0,27r] and [0"™,0f“] C [0,2?!]. 


3.1 





Figure 8. Given Oi, the LRL paths when the arc angie in the right turn is n. This figure shows the 
angies for 02 when the LRL path does not exist. 


Lemma4.1. := mm{ds,dR,RS^{ef"),RS^{df’^),RS^{e^^)}. 

Proof. Refer to the appendix for the proof. 

4.2. Optimizing the RL path 


□ 


We use similar notations as in previous subsections (refer to Fig. lOl. The angles 0 ( 6 i) and d2{RL, 0 i) 
are also written as (j) and 62, for brevity. The length of the RL path is denoted as T)(6i) and is equal to 
p{6i + 62 + 2^). RL paths do not exist when x > 4 p. In addition, even when 0 < x < 4 p, there are a subset 
of angles of di for which an RL path does not exist. Moreover, given di, there are two possible RL paths, 
as either <p + 62 < 7i or ^ + 62 > 7i. In the following discussion and in Fig. [T^ we assume that (l> + 62 < 7 i. 
The other RL path can be addressed similarly. 

We first define some values of 61 where the optimum can occur (these correspond to the extreme values 
of T) and 62 for the RL path and will be derived later in the proof). Let 6^* be the solution to the equation 
d2{RL, 61) = 01. Also, let 6^* and 6^* be the solutions to equation (/)(6i) + d2{RL, 61) = 71. Let dt denote 
the length of the shortest circular arc of type L that joins the two targets such that the boundary angles of the 
arc belong to the corresponding intervals at the targets and ^ + 62 < If such an arc does not exist, then di 
is set to 00. Let RL* = mm{RL^{d'f‘^^),RL\ef'‘),RL^{ef^'‘),RL^{df‘‘^}. 

Lemma 4.2. If x > Ip, minejg/j{RLi(ei)} := rmn{RL\d^*),RL\d^*),RL\d^*),RL*}. If0<x<2p, 
mine,g/j{RL^(0i)} :=rmn{dL,dR,RL^{e^*),RL*)}. 


Proof Refer to the appendix. 


□ 


5. Numerical results 


Computational results are presented for 25 instances with 20 targets in each instance. The locations of 
the targets were sampled from a 1000 x 1000 square. The minimum turning radius of the vehicle was chosen 
to be 100 . The heading angles at each target are discretized into 4 , 8, 16 , and 32 intervals. We use the Noon- 
Bean transformation to first convert the one-in-a-set TSP into an ATSR Then we use a transformation method 







outlined in to convert the ATSP into a symmetric TSP. This method converts an asymmetric instance 
with n nodes into a symmetric instance with 3n nodes. We chose this method primarily because unlike other 
transformations, there is no big-M constant involved, and therefore we did not have any numerical difficulties 
such as those faced in iT^T^^. For example, the transformed TSP instance with 32 discretizations 
at each target has 1920 nodes. Each of the transformed TSP instances was solved to optimality using 


the CONCORDE solver |24|. The improvement of the lower bounds as the number of discretizations or 


intervals increases is shown in Eig. 11 On average, the improvement of the lower bounds with respective to 


the optimal ETSP cost for 32 intervals was 22.28%. 

A feasible solution was also obtained by discretizing the angles at each target (32 values) and applying 
the above transformation procedure. The comparison of the cost of the feasible solution with respect to 
the optimal Euclidean TSP cost and the lower bound (corresponding to 32 intervals at each target) for 
the 25 instances is shown in Pig. The average deviation of the cost of the feasible solution from its 


corresponding lower bound is 5.2%, while the average deviation of the cost of the feasible solution from 
its corresponding ETSP cost is 29.2%. In one of the instances, we found the cost of the feasible solution 
from its corresponding lower bound improved by approximately 44%. These results show that the proposed 
approach can be used to obtain tight lower bounds for the DTSP. A feasible DTSP solution and an optimal 


solution corresponding to the lower bound for an instance are shown in Pig. 13 


6. Conclusion 


We provide a systematic procedure to find lower bounds for the DTSP. This article provides a new direc¬ 
tion for developing approximation algorithms for the DTSP. Currently, the transformation method increases 
the size of the one-in-a-set TSP by 2 or 3 times, resulting in a large TSP. Computationally, more efficient 
tools for directly solving the one-in-a-set TSP will be useful in finding fighfer lower and upper bounds for 
Ihe DTSP. Pufure work can also address fhe same problem wifh multiple vehicles and ofher precedence 
consfrainfs. 











Figure 10. RL path 



Instance Number 


Figure 11. Lower bounds computed with 4, 8,16, and 32 intervais at each target for 25 instances. 
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Figure 14. RS path: examples illustrating 2)(0i) and d 2 {RS,di). 






7. Appendix 


7.1. Proof of Lemma [id] 

Using Fig. [^ one can relate L and (/> to 6 i using the following equations: 

psin0 +Lcos(/) =^ —psinQi, 

pcos(/) — Lsin(/) = pcosFi. (20) 

The arrival angle 62 {RS, 6i) at target 2 is equal to27i — (p. We now consider two different cases: x>2p and 
X <2p (the RS path does not exist for a subset of angles of 6i if Jc < 2p). 

Case 1: X > 2p. 

The length of the RS path is D := (0i + (/))p +L. Therefore, ^ := (1 + ^)P + J^- The derivatives 
of 0 and L with respect to Q\ can be obtained by differentiating (|20l) as follows: 


, , T ■ j. 

(pcosp -Lsinp)-— +cosp-— = 
dd\ dd\ 

/ ■ JL r ^ dL 

-(psin^()+Lcos(/))—-sin(/)— = 


P cos 01, 
p sin 01 . 


Solving these equations and simplifying further, we obtain the following: 

^cos0 — 1, 
px 

— — cos 01 . 

Therefore, 


dcj) 

dL 

d6i 


dD d^ dL 

X 

= - (p cos (j) — P cos 01) 
= xsinip. 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 
(27) 


For any 0i G [0,27r], it is easy to verify geometrically that <j) G [0, tt] using Fig. [^ Therefore, V0i G (0,27r), 
^ > 0, i.e., the length of the RS path increases monotonically from Jc. When 0i = 271, the curved segment 
in the RS path vanishes and the length of the RS path returns to the Euclidean distance between the targets 
(x). Even though the length of the RS path increases monotonically for any 0i G [0,27r), the arrival angle 
at target 2, 02 := 271 — (j), first decreases with 0i, reaches a minimum at some 0i = 0*, and increases to 271. 
This minimum can be computed by solving = 0 |cos(0(0*)) — 1 = 0 or cos(0(0*)) = 4. One can 

verify that at 0i = 0*, 02 reaches a minimum. 

Now, the optimum for min^jg/j (0i)} must satisfy one of the following conditions: 

E g = 0 or 01 = 0 (g does not exist at this point) or 0i = 0for 0i = 0™“. V0i G (0,27r), g / 0. As 
the length of the RS path increases monotonically with respect to 0i, we need not consider 0i = 0™“^. 
Therefore, for this condition, the optimum occurs when 0i = 0 or 0i = 0J”'”. 


2. 02 = e"”'" or 02 = 

Therefore, whenT > 2 p, minejg/,{/? 5 '(ei)} := min{ds,IiS^ 


Case 2: X < 2 p. 

In this case, the RS path is not defined for any 0 \ G (sin(^),| +cos(^)). Moreover, when 0 \ = 
sin(^) or 01 = I +cos(^), the RS path reduces to just one segment of type R. Therefore, following the 
same analysis as in the previous case, mine|g/j{/? 5 '*( 0 i)} := min{(i5,dR,/?5'^(0™"),/?S^(0™”),7?S^(0““)}. 
Hence this case is proved. 

7.2. Proof of Lemma [ 4 ^ 

We can solve for (j) and 62 using the following equations (Fig. m- 


2 p cos (j) — p cos 02 = p cos 61, 

2psin0+psin02 =-^ —psin0i. ( 28 ) 

Differentiating and simplifying these equations, we get 


^ ^ dd) . ^ d02 . ^ 

-2smp-— + sme2^7r = -sm0i, 
d0\ d0\ 

^ ^ d(I) ^ d 02 

2 cos (p + cos 02 —r “ cos 01. 

d0i d0\ 


Solving further for and we get 




d(j) 

d 0 i 

d 02 

d 0 i 

dT) 

d 0 i 


sin(6i — 62) 
2 sm{(j) + 62) ’ 
sin(ei + 0) 
sin((/) + 62) ’ 


= p(l + 


d 02 

d 0 \ 


+ 2 


d^ 

d 0 \ 


sin(0i+0) sin(ei-02) s^ 
sin(0 + 02) sin(0 + 02) 


( 29 ) 

( 30 ) 

( 31 ) 

( 32 ) 

( 33 ) 


Equating = 0 and simplifying the equations, we get either 0 + 01 = 0 or 0 + 02 = 0 or 0 i = 02. 
0 + 01 = 0 or 0 + 02 = 0 would imply that one of the circles vanishes; however, this is possible only 
when Jc < 2 p. When 0 i = 02, we note that ^ = —1 and ^ = 0 . Using this, one can verify that 

~ ^ ^ Therefore, the length of the RL path reaches a minimum when 0 i = 02- 


Case 1: 4p > > 2p. 


The optimum for mine|g/j{/?L^( 0 i)} must occur at one of the extreme values of T)( 0 i) or when 0 i G 
{ 0 }""’, 0 j"“} or 02 G { 0 ™", 0 ™“}. T)( 0 i) reaches a local minimum at 0 i = 0 ^* (Fig. 151 . Also, the RL path 

or 01 = 0 ^*. Specifically, for a small £ > 0 , the RL path does not exist when 


-\min nmax 
^2 ’ ^2 

just ceases to exist when 0i = 0^* 








RL Length 




6i 

(a) x> 2p (b) jf > 2p 



Figure 15. RL Path: examples illustrating S)(0i) and 92{RL,ei) for the case when 0 < 0 + 02 < 


01 = 02* _£ or 01 = 02* + £. Therefore, min 0 , 6 /,{/?Li( 0 i)} := min{/?L^(0i*),7?L^(02*),/?L^(03*),7?L*}. 
Case 2: 2p > T > 0. 

In this case, one of the circles may cease to exist, and therefore the optimum may be equal to dt or dR 
if the corresponding angle constraints are met. Following the same arguments as in the previous case, we 
obtain mine|g/j{/?L^(0i)} := min{dL,dR,RL^ {6^*),RL*)}. Hence this case is proved. 













